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Abstract 

We evaluate the mass shift of isospin-averaged D-meson in the nuclear medium. 
Borel-transformed QCD sum rules are used to describe an interaction between the 
Z)-meson and a nucleon by taking into account all the lowest dimension-4 operators 
in the operator product expansion (OPE). We find at normal matter density the 
D-meson mass shift is about 10 times (~ 50 MeV) larger than that of J/ip. This 
originates from the fact that the dominant contribution in the OPE for the D-meson 
is the nucleon matrix element of m c qq, where m c is the charm-quark mass and q 
denotes light quarks. We also discuss that the mass shift of the -D-meson in nuclear 
matter may cause the level crossings of the charmonium states and the DD threshold. 
This suggests an additional mechanism of the J/ip suppression in high energy heavy- 
ion collisions. 
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1. Introduction 



Changes of hadron properties in the nuclear medium have recently attracted great interests 
in theoretical studies. They also have induced the on-going experiments and forthcoming 
experimental plans in heavy nuclei in GSI and in relativistic heavy-ion collisions in SPS at 
CERN and AGS, RHIC at BNL0. In particular the spectral changes of vector mesons are 
expected to be a possible signal to investigate such medium effect, because their leptonic 
decay in the medium can supply us information of the medium without disturbance of the 
strong interaction. Motivated by this experimental advantage, the in-medium effect of the 
light vector-mesons (p°,uj and 0) has been mainly studied in effective hadronic mo dels! 
and QCD sum rules (QSR's)IH§. Recently we applied the QSR method to low-lying heavy- 
quarkonium (J ftp) in order to investigate its mass shift in nuclear matter and the J/if>- 
nucleon interaction at low energy^. It was shown that the J /if) mass at normal matter 
density drops by 0.1 ~ 0.2% of its vacuum value. This negative mass-shift of J /if) corresponds 
to a small decrease of the mass, 4 ~ 7 MeV00. This is much smaller than the similar effect 
in the light mesons such as p° and to This difference can be understood as follows: In 
QSR the magnitude of the spectral change is related to the in-medium change of quark 
and gluon condensatesEUil. In the linear density approximation, the quark operators give 
a dominant contribution for p° and u, and have significant change in the matter. On the 
other hand, the dominant in-medium gluon condensates at the normal matter density for 
J /if) are only 5 ~ 10% smaller than its vacuum value. 

In this paper, we generalize the above calculations and apply the QSR analysis in Ref. |6| 
to light and heavy quark systems with unequal mass. We focus on in-medium properties of 
pseudoscalar mesons -D's. Studying the mass modification of the D-mesons in the nuclear 
medium is important by the following physical reasons: 

1. A part of total cross section of the D-meson production attributes a conversion of 
the cc-states such as J /if) and if)' into open-charm pairs. The component is produced 
through the reaction, if) + N ^ D + D, in high energy proton- nucleus collisions and 



relativistic high energy heavy-ion collisions. It contains a heavy charm-quark and a 
light quark. The existence of a light quark in the D-meson causes much difference of 
the properties in the nuclear medium between the Z)-meson and J/ip. The latter is 
dominated by the gluon condensates as discussed above, while the former has a large 
contribution from the light-quark condensates, multiplied by a charm-quark mass. 
Furthermore J /if) predominantly interacts with light hadrons solely through gluonic 
content effects in matter, while the D-mesons will couple strongly through inelastic 
channels such as DN — > A c or S c (+7r). Thus we can expect a large modification of 
the D-meson in the nuclear medium as well as p° and u. 

2. The existence of novel states, so-called the D-mesic nuclei, has been predicted by 
quark meson coupling (QMC) modelB The QMC model suggests that D~(cd) may 
be bound in heavy atoms such as 208 Pb by an attractive scalar-meson exchange and an 
attractive Coulomb force. It also suggests that D°(cu) is deeply bound in the nuclei 
by an attractive tu-meson exchange. In view of QSR, it is of importance to investigate 
the D meson-nucleon (N) interaction. 

3. The nuclear absorption of the D-mesons created in ir-A and p-A collisions (A denotes 
targets such as Be, Cu, Al, W and Au) has been measured via the decay D — > Kir 
at Fermilablll. The result indicates that D-mesons are not completely absorbed in the 
nuclear targets irrespective of the charge of the produced .D-mesons. The theoretical 
analysis for the D-N interaction will give an important suggestion for the experimental 
data in this case too. 

Motivated by these points, we investigate the properties of .D-mesons in the nuclear matter 
through the D-N interactions, using an application of Borel-transformed QSR to the D- 
N forward scattering amplitude. Here we deal with isospin-averaged D-meson current for 
simplicity. 

After a brief explanation of the QSR formulation, in the next section, we calculate the 
mass shift at finite density using the Borel QSR and compare the results with the mass 
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shift of J/i[) calculated in Ref. [10]. On the basis of these results, we will discuss a possible 
mechanism of J/ip suppression through the level crossings between the charmonium states 
and the DD threshold. 



2. QCD sum rule analysis for D-meson mass shift 

We start with a two-point in-medium correlation function Hp^g to discuss hadron properties 
in the nuclear matter. In the Fermi gas approximation for the matter, Ilp^f is divided into 
two parts by applying the operator product expansion (OPE) to the correlators in the deep 
Euclidean region (Q 2 = —q 2 > 0). One is a vacuum part, Upg, and another is a static one- 
nucleon part, Tp$- This decomposition is expected to be valid at relatively low density^. 
Namely, in the framework of QSR, IIp^ can be approximated reasonably well in the linear 
density of the nuclear matter that all nucleons are at restHPHS: 

O?) = i\<£x e i ^(TJ 5 (a;)4(0))NM(p iV ) * H° Pg (q) + ^T PS (q), (1) 

where pN denotes the nuclear matter density, and the forward scattering amplitude Tps of 
the pseudoscalar current-nucleon is 

T PS (u;,q) =i J d 4 xe^ x (N(p)\TJ 5 (x) 4(0)\N(p)). (2) 

Here = (ui,q) is the four-momentum carried by the .D-meson current, J${x) = j\{x) 
= (cij5q(x) + 575c(rr))/2, where q denotes u or d quark. \N(p)) represents the isospin, 
spin-averaged static nucleon-state with the four-momentum p = (M^,p = 0), where 
is the nucleon mass, 0.94 GeV. The state is normalized covariantly as (N(p)\N(p')) = 
(27i) 3 2p° 5 3 (p — p'). The QSR analysis on the forward scattering amplitude enables us to 
obtain the information for the D-N interaction. In Eq. ([]]), the second term means a slight 
deviation from in- vacuum properties of the D-meson determined by n P5 . By applying QSR 
to Tps-, we get the D-N scattering length ap> in the limit of q — > 0. In this limit, Tp$ can be 
related to the T-matrix, T v x{m D) q = 0) = 8ti(M n + m D )aD- Near the pole position of the 
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D-meson, the spectral function p(u, q = 0) is given with three unknown phenomenological 
parameters a, b and c in terms of the T-matrix: 

%>tf(w, 0) 



p(-,q = 0) = -i^lm 

it <m £ 



(uj 2 — m 2 n + is) 2 



+ ■■ (3) 



D 

= ad'(uj 2 - m 2 D ) + b5(uj 2 - m 2 D ) + c5(uj 2 - s ). (4) 

Here the leptonic decay constant fo is defined by the relation (0\J 5 \D(k)) = foTn 2 D /m C) 
where \D(k)) is the D-meson state with the four-momentum k. rrip, = 1.87 GeV and 
m c = 1.35 GeV are masses of the D-meson and charm quark respectively. The terms denoted 
by ■ ■ • in Eq. (|3|) represent the continuum contribution and 5' in Eq. (^) is the first derivative 
of 5 function with respect to uj 2 . The first term proportional to a is the double-pole term 
corresponding to the on-shell effect of the T-matrix and a is related to the scattering length 
cld as a = —8n(MM + rnD)aDfD m D/ m c- The second term proportional to b is the single-pole 
term corresponding to the off-shell effect of the T-matrix. The third term proportional to c 
is the continuum term corresponding to other remaining effects, where So is the continuum 
threshold in the vacuum. Combining the single-pole term of Ilp 5 in Eq. (|I|) with Eq. 
we can relate the scattering length extracted from the QSR of Tps with the mass shift of 
the .D-meson, 

M N + m D 

dm D = 2tt— p N a D - (5) 

M N m D 

We may determine these unknown parameters by using a dispersion relation and matching 
the phenomenological (ph) side with the OPE side. Before such analysis, we can impose a 
constraint from the low energy theorem among these parameters: In the low energy limit 
(w -> 0), T ps {uj,0) becomes equivalent to the Born term T|° rn (cj,0), Tj£(0) = T|° rn (0). 
We assume two cases to take the Born term into the "ph" side. The case (i) and the case (ii) 
are defined as followJl. At g M ^ 0, we require the a; 2 -dependence of the Born term explicitly 
in the case (i): 

Tf s {u 2 ) = T|° m (^ 2 ) + , 2 a 2 , 2 + + — ^-S. (6) 

(mf) — u 2 ) 2 m z D — uj 2 Sq — ui 2 



with the condition 

a b c /«\ 
— + — + - = 0. 7 

In the case (ii), we do not require the a; 2 -dependence of the Born term explicitly: 

T ^ s{uj2) = (m* D -u*)* + m* D -u* + 7^?' (8) 
with the condition 

^ + ^ + f = ^ B no). o) 

m D m D S 

If Tp™ n (0) = 0, two cases coincide. We determine two unknown phenomenological param- 
eters a and b in the QSR after the parameter c is removed by the condition of Eq. ([?D or 
(§). The analysis of the Born term is easily performed through a calculation of the Born 
diagrams at the tree level. The isospin states of the D-meson determine the contribution 
to the Born term as follows: For the D°(cu)-N and D + (cd)-N interactions, we need two 
reactions, 

D°(cu) + p(uud) or n(udd) — ► A+, £+(cwd) or E° c (cdd) (10) 
and 

D + (cd) +p(uud) or n{udd) — ► £+ + (cuu) or A+,£+(cu<i), (11) 



where the static masses of charmed baryons are 2.28 GeV for A+ and 2.45 GeV for S°, S c 
and Note that we use M B ~ 2.4 GeV as the average value, where B means either A+, 

E+, £+ + or E°. After averaging over the nucleon spin, we obtain 



T->Born / Q) _ 1 M N (M N + M B ) 

The pseudoscalar form-factor g p (q 2 ) is introduced through the following relation^, 

(B(p')\Qi l5 q\N(p)) = g P (q 2 Mp')z l5 u(p), (13) 
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where u(p) is a Dirac spinor and Q(q) in the pseudoscalar current is the heavy(light)-quark. 
Furthermore, g p {q 2 ) is related to a coupling constant gNDB such as 

2 2m 2 D f D g NDB , A 

0p(? = ; 2 2"' 14 

where m q denotes a light-quark mass. gND\ c has been estimated in Ref. O which gives 
gDNA c — 6.74. However since gDNT, c has not been evaluated, we take an approximation 
gDNH c — gNDk c - On the other hand, there are no inelastic channels like the above for 
D°(cu)-N and D~(cd)-N interactions, i.e. Tj^ rn (0) = in both case (i) and case (ii). 

Eventually we determine parameters a and b simultaneously after applying the Borel 
transform to both the OPE side and the "ph" side. 



3. Numerical results in the Borel sum rule 

After performing the Borel transform B of the Tp$ calculated in the OPE up to dimension- 
, we obtain as a function of the Borel mass M 



B 



T, 



OPE 
PS 



I 777, 

-m c (qq) N + 2{qHD q) N ( -1 + j± 



o m c 
M 2 



+ as ((v r\ 2 1 r 2 \ 1 1 4 lm l + l ( m V 



+e m 2 JMl I _ 



""C -1 
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(15) 



where (-)jv denotes the nucleon matrix element. The renormarization scale is taken to be 
fi 2 = 1 GeV 2 and the Euler constant is 7^ = 0.5772 • ■ ■. Corresponding formula for the "ph" 
side, in the case fi) reads 



B 



T, 



ph 



PS 



a f_L e X/^ - +6 | ( -r, 

\ M 2 171% J 



m 2 D / M 2 _ S c _ sq/M 2 

m 2 D 
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1 M N (M N + M B ) 4f 2 D m% 2 
2(M N + M B ) 2 - m 2 D (m q + m c ) 2 9ndb 



x 



' 1 1 —m D /M 



{M N + M B ) 2 - ml \{M N + M B ) 2 - ml M 2 
and that in the case (ii) reads 



e 



(16) 



B 



T ph 

J PS 



\M 2 m% J \ m D 



+s T|^ rn (0)e- so/M2 . (17) 



In Eq. (|15D , a convention, (it • G) 2 = G" A Gp A m k m p is used, where u = (1, 0) for the static 
nucleon. We equate Eq. fllSD with Eq. (|16D or Eq. (|I7D in the case (i) or the case (ii) 
respectively. Furthermore we take a first derivative of its equation with respect to M 2 in 
each case. We can derive the Z)-meson mass shift as a function of M 2 by removing the 
parameter b from two equations obtained thus. The Borel curve for the D-meson mass shift, 
Smo, is shown in Fig.l. Here we adopt fjj ~ 0.18 GeV for the decay constant and so — 6.0 
GeV 2 for the continuum threshold. The values of these parameters were estimated by the 
analysis of the in- vacuum correlation function in Ref. [n|^U[ This fu value is very close to 
the result of the lattice QCD calculation (0.194 ± 0.01 GeV)0 and is consistent with the 
upper bound of experimental data (< 0.31 GeV)il. Note that the parameters are fixed in 
our calculation. We use the nucleon matrix elements^! for quark fields such as {qq)N — 5.3 
GeV and (qUD q) N ~ 0.34 GeV 2 , and for gluon fields such as (^G 2 ) N ~ -1.2 GeV 2 and 
^((u ■ G) 2 - \G 2 ) N ~ -0.1 GeV 2 . The light-quark mass is taken to be m q ~ 0.008 GeV. In 
Fig.l, the solid line ("Born = 0" [case (i)]) are calculated without the contribution of inelastic 
channels. In the doted line ('With Born term" [case (ii)]) and the dashed line ("Without 
Born term"), we allow for the contribution of the sub-threshold resonances A c and S c , lying 
very close to the D-N threshold. As shown in Fig.l, the sub-threshold effect is repulsive 
both in the case (i) and the case (ii) and its contribution in the case (ii) is stronger than 
that in the case (i). This is analogous to the case for the K~-p inter act ioniH. In Fig.l, 
the remaining three lines are evaluated by extracting the contribution of m c (qq)N term in 
the OPE from the solid, dashed and doted lines respectively. We find that the contribution 



of m c {qq)N term is more than 95% of total contribution corresponding to the solid, dashed 
and doted lines within the plateau regions. 

The analysis of this graph is summarized in Table 1. We cannot determine th Borel 
windows, since there is only lowest-dimension term in the OPE side. Therefore we take the 
following procedure to determine a window of M 2 : We focus on a plateau region of each 
line shown in Fig.l. First, we take the minimum point in the line as the smallest value of 
Snip- Next, we determine two points in M 2 as the deviation from the minimum value of 
Stud becomes less than 10% of the minimum value. We take the region between the two 
points as a window of M 2 . The window of M 2 in each line is given in Table 1. As shown in 
Table 1, all the Borel curves are rather stable within the windows. Note that the windows 
discussed above are also close to that obtained by scaling up typical Borel-windows for the 
light-vector mesonsi. An estimate from the light-vector mesons is (^-) 2 x 0.8 (= 4.7) < 
M 2 < (^) 2 x 1.3 (= 7.7) for V = p, uj and (^) 2 x 1.3 (= 4.4) < M 2 < (^) 2 x 1.8 (= 6.1) 
for V = 4>. 

It should be stressed that the QSR for vacuum correlation function H° PS up to dimension- 
6 operators cannot reproduce well the D-meson mass (mo) in free spaced. The Borel curve 
of mo at So = 6.0 GeV 2 does not have any stability and seems to give larger values than 
the experimental data (m_o = 1.87 GeV) in the plateau region (M 2 = 3 ~ 8 GeV 2 ) of 
5rri£) discussed above. Furthermore, the curve of has rather larger change than that 
of 8m d in the region. Therefore, if we perform the QSR analysis for the effective mass 
m* D = vhy) + SrriD, the stability in the Borel curve of m* D will become obscure. This implies 
that the QSR analysis for m* D is not so valid in this case. However the application of the 
QSR to SrriD gives us a possible indication for the D-meson mass-shift as discussed above. 

From the above analysis we obtain 5m d = —48 ± 8 MeV by considering both cases ((i) 
and (ii)). Then the D-N scattering length is —0.72 ± 0.12 fm. This result suggests that 
the D-N interaction is more attractive than the J/ip-N interaction, where Sm j/^ is about 
— 5 MeV and aj/^ is about —0.1 fm00. It is noted the above results were performed at the 
nuclear matter density, = p ~ 2p , where p is the normal matter density. Assuming 
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that the linear density approximation in the QSR is valid at the nuclear matter density, we 
find these results lead to a larger decrease of the D-meson mass than that of the charmonium. 
We also find that the large mass-shift of the D-meson originates from the contribution of 
m c (qq) n term in the OPE. In the QMC model, the contribution from the m c (qq) n term may 
correspond to a quark-cr meson coupling. In fact, the model predicts the mass shift of the 
D-meson becomes —60 MeV for the scalar potential at the normal matter density^. Their 
results are very close to our results. 

4. Concluding remarks 

We present an analysis of isospin-averaged Z)-meson mass-shift through a direct application 
of the Borel QSR to the forward scattering amplitude of the pseudoscalar current and nu- 
cleon. Here we perform the operator product expansion with all the terms up to dimension-4. 
The result predicts an attractive mass-shift of the isospin-averaged D-meson about 50 MeV 
(about 3% of the bare mass) at the normal matter density. The contribution of m c (qq) jy term 
in the OPE is more than 95% of the results evaluated up to all the dimension-4 operators. 
Our result is very close to that reported by the QMC mo delR Then the D-N scattering; 
length do is about —0.7 fm and indicates a strongly attractive force between the D-meson 
and nucleons. The mass modification of the D-meson in the nuclear matter corresponds to 
10 times larger than that of J/if>. 

Recent NA50 Collaboration^ has reported a strong (so called "anomalous") J /if) 
suppression^ in Pb-Pb collision at 158 GeV per nucleon. We suggest a following possi- 
bility to cause such a strong J/ if) suppression using the above large difference of the mass 
shifts between the D-meson and J /if): Suppose that the other cc-states such as if)' and Xc's 
do not also have a large mass-shift in the nuclear matter as discussed above for J /if). Then 
on the basis of our calculations, the DD threshold (~ 3.74 GeV in free space) decreases by 
100 MeV at p^ = Po and comes down between Xc2 (~ 3.55 GeV) and if)' (~ 3.68 GeV). 
At p N = 2p , it decreases by 200 MeV and comes down between Xci (~ 3.50 GeV) and 
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Xc2 (~ 3.55 GeV). Such level crossings between the DD threshold and the charmonium 
spectrum are shown in Fig.2. It is well known from the p-A collision data that only 60% of 
J /if) observed are directly produced and the remainder comes from excited states (% C (1P), 
■^'(25')) with a ratio of 3 to 10. So if the disappearance of if)' and Xc due to their decays 
into DD takes place, it could lead to a decrease of the J /if) yield. This J /if) suppression 
will have stair-shaped form as a function of the energy density of the system. At low energy 
density, only if)' — > DD takes place, which causes a slight suppression of J/ip. At interme- 
diate density, subsequent suppression occurs by the level crossing of (xc2, Xa) states with 
the DD threshold. At higher density, direct suppression of J/ip by the decay J /if} — ► DD 
occurs. This is an alternative or supplementary effect to the deconfmement scenario in Ref. 
29] and may have implications to the recent data of NA50 Collaboration!!. Strictly speak- 



ing, as mentioned before, the formulation of QSR used here is applicable to dilute nuclear 
matter in comparison to highly dense nuclear matter produced in high energy heavy-ion 
collision. If we try to discuss the mass shifts at higher density, it will need an extension or 
some improvements in the above formulation. This is one of important future works for us. 
However we may say that our calculation in this paper can suggest a possibility of the level 
crossings even in only 1 or 2 times normal matter density and lets one expect the direct J / if) 
suppression at sufficiently higher density. 

We may also suggest another phenomenon caused by the level crossings discussed above. 
It is a change of the decay width of the charmonium. In the vacuum the resonances above 
DD threshold, for example the if)" state, have a width of order MeV because of the strong 
open-charm channel. On the other hand, the resonances below the threshold have a very 
sharp width of a few hundreds keV. So after the level crossings, the decay modes of the if)' 
and Xc states will change drastically at least one order of magnitude. 

Needless to say, in the high energy heavy-ion collisions we must also perform the theo- 
retical investigation of finite temperature effecJE'H to the mass modification. Therefore, 
a future task will be to understand the medium effect in terms of the QSR when both 
temperature and density are finite. As for investigation of only the D-N interaction, an 
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inverse kinematics experiment will give useful informational Since the projectile (target) is 
a heavy-ion beam (light-nuclei), the decays of the D-mesonEl and the charmonium inside a 
nucleus will be possible in such an experiment. 

Before closing, it is stressed that we should take the above analysis for the isospin- 
averaged D-N interaction as a qualitative estimate. If we wish more quantitative discussion, 
the isospin-decomposition of the D-meson and higher correction terms beyond dimension-4 
operators must be taken into account. Particularly the odd-components in the OPE play 
important roles for a difference of the mass shift between D and D. 
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FIGURES 

Fig. 1. The L>-meson mass-shift 5m d [GeV] at normal matter density (po = 0.17 fm -3 ) as a 
function of the squared Borel mass M 2 [GeV 2 ] : The solid, doted and dashed lines are calculated 
by taking into account all the dimension-4 operators in the OPE. The lines correspond to the cases 
of 7fg rn (0) = 0, with (case (i)) and without (case (ii)) an explicit Born term (T^ rn (0) + 0) in 
the "ph"side respectively. The remaining three lines show the results obtained by extracting the 
contribution of m c (qq)N term in the OPE from the above each line respectively. 

Fig. 2. Comparison between charmonim spectrum and DD threshold level: If the mass shifts of 
all the charmonium states are very small (0.1 ~ 0.2% of the bare mass) in the nuclear medium 
(pN = Po ~ 2po ; Po = 0.17 fm~ 3 ), we predict that the threshold lying just below ID state in free 
space falls down below the ip' state at normal matter density and the Xc2 state twice at the density. 
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TABLES 

Table 1. The result of an analysis for the solid, dashed and doted lines in Fig.l: The second raw 
shows the plateau region, which we determine as the deviation from a minimum value of Smo in 
each line becomes less than 10% of the minimum value. The 5m d in the plateau region is given in 
the third raw. 



Lines 


Solid line 


Doted line (case (i)) 


Dashed line (case (ii)) 


Plateau region [GeV 2 ] 


3.5 < M 2 < 7.0 


2.6 < M 2 < 7.0 


3.4 < M 2 < 8.0 


6m D [GeV] 


-63 ±3 


-53 ±3 


-44 ±3 
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